Python에서 직접 Simulink 모델 시뮬레이션하기

function res = sim_the_model(args)
% 지정된 매개변수로 Simulink 모델을 시뮬레이션하는 유틸리티 함수
% 
% 입력:
%    StopTime: 시뮬레이션 종료 시간, 기본값은 nan
%    TunableParameters:
%       필드가 조정 가능한 참조 워크스페이스 변수이며 시뮬레이션에 사용할 값을 가진 구조체.
%    
%    위 입력에 대해 nan 또는 빈 값은 시뮬레이션이 모델에 설정된 기본값으로 실행되어야 함을 나타냅니다.
% 
% 출력:
%    res: 기록된 신호의 시간 및 데이터 값을 가진 구조체.

arguments
    args.StopTime (1,1) double = nan
    args.TunableParameters = []
end

    %% SimulationInput 객체 생성
    si = Simulink.SimulationInput('suspension_3dof');
    %% SimulationInput 객체에 StopTime 로드
    if ~isnan(args.StopTime)
        si = si.setModelParameter('StopTime', num2str(args.StopTime));
    end
    
    %% 지정된 조정 가능한 매개변수를 시뮬레이션 입력 객체에 로드
    if isstruct(args.TunableParameters) 
        tpNames = fieldnames(args.TunableParameters);
        for itp = 1:numel(tpNames)
            tpn = tpNames{itp};
            tpv = args.TunableParameters.(tpn);
            si = si.setVariable(tpn, tpv);
        end
    end

    %% sim 호출
    so = sim(si);
    
    %% 시뮬레이션 결과 추출
    % 기록된 신호의 시간 및 데이터 값을 구조체에 패키징
    res = extractResults(so,nan);

end % sim_the_model_using_matlab_runtime

function res = extractResults(so, prevSimTime)
    % 기록된 신호의 시간 및 데이터 값을 구조체에 패키징
    ts = simulink.compiler.internal.extractTimeseriesFromDataset(so.logsout);
    for its=1:numel(ts)
        if isfinite(prevSimTime)
            idx = find(ts{its}.Time > prevSimTime);
            res.(ts{its}.Name).Time = ts{its}.Time(idx);
            res.(ts{its}.Name).Data = ts{its}.Data(idx);
        else
            res.(ts{its}.Name).Time = ts{its}.Time;
            res.(ts{its}.Name).Data = ts{its}.Data;
        end
    end
end

function figHndl = plot_results(res, plotTitle)
%PLOT_RESULTS call_sim_the_model의 결과를 플로팅

figHndl = figure; hold on; cols = colororder;


plot(res{1}.vertical_disp.Time, res{1}.vertical_disp.Data, 'Color', cols(1,:), ...
    'DisplayName', '수직 변위: 기본 Mb 값으로 첫 번째 시뮬레이션');
plot(res{2}.vertical_disp.Time, res{2}.vertical_disp.Data, 'Color', cols(2,:), ...
    'DisplayName', '수직 변위: 새로운 Mb 값으로 두 번째 시뮬레이션');


hold off; grid; 

title(plotTitle,'Interpreter','none');
set(get(gca,'Children'),'LineWidth',2);
legend('Location','southeast');

end

si = simulink.compiler.configureForDeployment(si);

appFile = which('sim_the_model2');
outDir = fullfile(origDir,'sim_the_model2_python_package_installer');
compiler.build.pythonPackage(appFile, ...
        'PackageName','sim_the_model2', ...
        'OutputDir',outDir);